Mean field superconductivity approach in two dimensions: Hydrogen in Graphite 
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Within the BCS theory of superconductivity we calculate the superconducting gap at zero tem- 
perature for metallic hydrogen-graphene system in order to estimate the superconducting critical 
temperature of quasi two dimensional highly oriented pyrolytic graphite. The obtained results are 
given as a function of the hydrogen-induced density of carriers n and their effective mass m*. The 
obtained gap shows a Maxwell-like distribution with a maximum of ~ 60 K at n ~ 3 x 10 14 cm -2 and 
m* /m = 1. The theoretical results are discussed taking into account recent experimental evidence 
for granular superconductivity in graphite. 
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I. INTRODUCTION 

Since the discovery of high temperature 
superconductivity we learned that low dimensional 
structures are good candidates for high critical supercon- 
ducting temperature T^. Related to this dimensionality 
effect it appears appropriate to mention also the role 
played by two dimensional (2D) interfaces in triggering 
superconductivity in nominally non-superconducting 
environment 3 -^ 5 -. One of the main paradigms for 2D 
carrier systems, graphene as well as highly oriented 
pyrolytic graphite (HOPG) or multigraphene have 
attracted considerable attention recently, partially due 
to the possibility of regulating the carrier density n 
by the field effect^. On the other hand, hints for 
the existence of superconductivity in HOPG have been 
published in the last years 8 ^. 

Let us first point out some interesting aspects of multi- 
graphene that should play a role in triggering the su- 
perconductivity phenomenon. Whereas graphene should 
be nominally a 2D system, multigraphene or HOPG 
is a quasi-2D system due to the weak coupling be- 
tween graphene layers, which is the reason for the huge 
anisotropy in resistivity, for example 10 . Both materials, 
graphene as well as HOPG have high energetic phonons 
going up to ~ 0.15 eV with many other branches at lower 
energie o n i 12 . High energy phonons due to small element 
mass provide an excellent condition for superconductiv- 
ity as in the case of metallic H at high pressures^ 3 .. From 
this point of view it seems reasonable to look for super- 
conductivity in these materials taking into account the 
whole phonon spectral function weighted by the electron- 
phonon coupling^. 

For HOPG it has been recently shown that 
upon temperature (and defect concentration) n ~ 
10 8 . . . 10 12 cm -2 J£ The lowest value obtained in clean 
samples and at low temperatures is clearly smaller than 
in free-standing as well as fixed-on-substrate graphene 
samples. To achieve a clear increase in the carrier den- 
sity the case of attached hydrogen in HOPG is of spe- 



cial interest. According to theoretical results upon the 
amount and the way hydrogen is attached to a graphene 
layer, a semiconducting, metallic or even a magnetically 
orderer state appear o 16 i 17 ' 18 . In case of metallization the 
electron density may increase dramatically in the region 
near hydrogen, e.g. 0.1 to 1 electron per unit cell. Added 
to this effect, the effective mass increases from the usual 
very small values m* < 0.01m to nearly the free electron 
mass nv*k Note that in this case no Dirac dispersion rela- 
tion is valid for graphite carriers but the usual quadratic 
one, dispersion that we will assume through all this work. 
The hydrogen-graphene bound system is interesting due 
to the influence in the spectral density similarly to the 
case of superconductivity in aromatic molecule a 14 ' 19 . 

On the other hand, we know that the reaction of hy- 
drogen adsorbed on carbon is endothermic and difficult 
to realize. However, experimental results indicate large 
amounts of hydrogen is present in HOPG samples 20 . Hy- 
drogen in HOPG might be included through the synthesis 
of this material, obtained after high-pressure and high- 
temperature treatment of polymers (e.g. kapton) foils. 
We expect therefore hydrogen can be bounded to carbon, 
specially around vacancies and other defects, or at the in- 
terfaces conforming regions or pockets with different elec- 
tron densities. Hydrogen bounded at defects might be the 
origin for the ferromagnetic propertie s 21 ' 22 , the anoma- 
lous transport observed in HOPG in the last yearaiS 
and for the existence of non-homogeneous, granular-like 
superconductivity- 9 , possible at the interfaces between 
highly crystalline graphite regions in HOPG samples 23 . 

The purpose of this work is to calculate the supercon- 
ducting energy gap at T = K for a graphene layer and 
extend this result to the anisotropic graphite case. Fol- 
lowing the Mermin- Wagner theorem in a pure 2D sys- 
tem there is no superconductivity^ and therefore the 
problem should be treated as a Beresinskii-Kosterlitz- 
Thoulcss (BKT) phase transition 25,2 ^ 5 . In the anisotropic 
case the problem is solved having a small anisotropy e, 
which resembles the low coupling between graphene lay- 
ers in graphite. As in the 2D anisotropic Heisenberg 
model, no matter how small is e, one has an appreciable 
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critical temperature because of the logarithmic behav- 
ior of the thermal fluctuation influcnc o 27 ' 28 . The results 
of this work indicate that hydrogen could play a deci- 
sive role triggering high temperature superconductivity 
in graphite. 



II. SUPERCONDUCTIVITY AT HIGH 
METALLIC DENSITY OF CARRIERS 

The estimate of the critical temperature T c using BCS 
is done from the energy gap equation given b y 29 i 30 



A(E) 



-N(0) J dE'V(E- 



E')A(E')(l-2f(E'))/2E', 
(1) 

where N(0) is the density of states at the Fermi energy 
E F , V(E - E') = V P {E - E') + V c is the interaction po- 
tential that we split into the electron-phonon term and 
the Coulomb potential, this last taken as a constant for a 
given carrier density. The potential Vp(E — E') depends 
on the pair-interaction- assisted phonon energy E — E' . 
Both potentials will be estimated below. Finally the 
Fermi-Dirac distribution function f(E') at a given T that 
we set at T c . 

At high metallic densities Ep 3> Ep> (this last the De- 
bye energy) T c is estimated assuming that electron pairs 
are formed with an energy difference up to Ep> around 
Ep. Then, the well known equation for the zero temper- 
ature energy gap (weak coupling limit) 14 i 29 i 30 

A(T = 0)«2£^exp(-l/(A-/x*)), (2) 

is obtained with A = N(0) < Vp > and 

M* =fi/{l+nln(E F /E D )), (3) 

with fi — N(0) < Vc > where the < ... > means the 
average within Ep or Ed and the screening length of the 
corresponding potentials. The well known result ^ is 
obtained solving Eq. (fT]) by introducing cutoffs energies 
for electrons and phonons. Using an Einstein approxi- 
mation Morel and Anderson obtained the same results 
for A(0),A and fi* without introducing cutoffs 3 -!. The 
results remain the same for a 2D system but due to the 
BKT phase transition it is valid only at T — K. 



III. ESTIMATE OF A(0) IN 3D AND 2D AS A 
FUNCTION OF CARRIER DENSITY n 



We discuss now the differences between 2D and 3D. 
We express all the necessary quantities as a function 
of r s , the distance between the carriers in units of the 
Bohr radius. The quantities qp(r s ) — 1.91/(ob'"s) and 
1t{t s ) — l-38/(a_B V /r7) are the 3D values for the Fermi 
vector and the inverse of the Thomas-Fermi decay length. 
Analogously for 2D we have {r s ap) and 2/ob. No- 
tice that the inverse decay length in 2D does not depend 



on the density n — (tt) 1 (asr s ) 2 . Another important 
difference is in the Coulomb interaction that behaves as 



for 3D and 
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for 2D, where q is the wave vector. These potentials 
should be reliable for r s < 15 above which Pines^ 2 . noted 
that convergence problems in the estimates may occur. 
Note that the value r s = 15 is still smaller than those 
needed for the Wigner crystal formatio n 33 i 34 . 

In the 3D case the values of A and /i in Eqs. @ and 
© are given by 
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r 2 ln((l + r 2 )/r 2 ) 



(6) 
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where r 2 = (qT(r s )/2qp(r s )) 2 . The Eqs. (6,7) corre- 
spond exactly to Eqs. (46,44) from Ref. 31. Note that our 
Eq. (7) is the average angular value of Eq. (26) of Ref. [3l| 
that missed a factor r 2 in the corresponding Eq. (44). 
For the case of 2D we have ii as the average angular of 
the Coulomb potential 
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and A is calculated as the average angular value of the 
square Coulomb potential multiplied by 0.66 because of 
averaging the q— moment. Therefore we have 



A = 0.66 
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Now we are prepared to estimate the values of A(0) in 2D 
and 3D for different values of r s , a function of the elec- 
tronic density. It should be notice that A and \i obtained 
as defined previously agree with the values obtained by 
Morel and Anderson^. They obtained them including 
an Einstein effective frequency as approximation for the 
phonon structure. Within a Debye model similar equa- 
tions are obtained but the phonon contribution is taken 
as an average on the whole phonon spectra. 

We found that A(0) is in general about 10 . . . 100 times 
larger in 2D than in 3D at similar values of r s . In Fig. 1 
we present A(0) vs. r s for the 3D case for 3 values of m* 
corresponding to 0.1, 1 and 3 the free electron mass. Fig- 
ure 2 shows the results for the 2D case. What can be the 
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FIG. 1: Energy gap at T = K vs. the distance between 
carriers r s in units of Bohr radius as = ATveoh 2 /rn* e 2 for 
the 3D case and for the effective masses m* — 3,1,0.1 of 
the free electron mass (continuous, dashed and dotted lines, 
respectively). The curves were obtained using Eqs. ([2]3,4,6,7) 
and the parameter Ed = &b860 [K] as an average over all 
frequencies. 



physical reasons for having energy gaps much larger in 
2D than in 3D? They are related to the different weights 
and values of the ratios qp/^qp and qr/qD (<Zr> is the 
averaged maximum phonon wave vector). Their contri- 
butions in 2D are dominant and depend on the effective 
mass, as seen in Fig. 3 where we plot the values of qr/^qp 
and qx/qD, with qn = 10 8 cm -1 for 2D (as in 3D). Note 
that the differences in A(0) between the 2D and 3D cases 
are larger the smaller the effective mass m* . 

It is adequate to study quasi 2D problems with a degree 
of anisotropy provided by the ratio between the conduc- 
tance in the atomic lattice plane and perpendicular to it, 
i.e. e ~ yo^Joab. For example, in the case of HOPG 
we have an in-plane conductance > 10 4 than the con- 
ductance perpendicular to the plane. As in the case of 
the Heisenberg model in 2D with anisotropy, even if it 
is small, this is very effective to recover again the super- 
conducting phase transition at T > K by a logarith- 
mic function of this anisotropy, i.e. T c oc l/ln(if(l — 
e)) oc — l/ln(e) for e — + 0, where K(x) is the elliptic 
functio n 27 ! 28 . The anisotropy produces a variation of T c 
of ~ 10 — 20% of the value obtained from the simple esti- 
mate using the BCS result T c ~ 2A(0)/(3.5fc B ). There- 
fore the value of A(0) remains a good approximation for 
the anisotropic quasi 2D case. 

We performed first calculations for HOPG at relatively 
low densities. We note that in graphene as well as in 
HOPG the expected T c is negligible small at electron 
densities n < 10 13 cm -2 . According to the BCS ap- 
proach and our estimates neither materials can be su- 
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FIG. 2: The same as in Fig. 1 but for the 2D case with the 
same parameters. 
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FIG. 3: Thomas-Fermi wavevector ratios qr/2qF = r s /V2 
(dashed-dotted line) and qr/qo for the 2D case for the three 
different effective masses 3, 1 and 0.1 of the free electron mass 
(continuous, dashed and dotted lines). 



perconductors at those electron densities. However, for 
the electron densities obtained for the case of hydro- 
gen fixed in special (e.g. Stone- Wales) defects we have 
n ~ 10 14 ... 3 x 10 14 cm~ 2 with a metallic-like band 
character that crosses the Fermi level with m* ~ mi 7 -. 
The estimated critical temperatures for these densities 
are around 25 K, values comparable to the observations 
in multigraphene samples^ 3 .. Note that decreasing n by 
a factor of two increases T c approximately by a similar 
factor, see Fig. 2. According to our estimate, a further 
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decrease of n, however, would produce a decrease of the 
critical temperature, see Fig. 2. 

Although defects in the graphene/graphite structure 
would increase the carrier density in a first stage, we be- 
lieve that attached hydrogen may be the most probable 
reason for the local increase of n. The idea that hydro- 
gen may trigger superconductivity in graphite due to the 
large increase in the electron density provides also a way 
to understand several experimental facts that we sum- 
marize below. Experimental results from SQUID£ ' 35 i 36 
and magnetotransport 9 -^ indicate that the possible su- 
perconducting state in HOPG has granular character: 
neither percolation in transport nor Meissner effect are 
observed so far for samples larger than ~ 10 /im. Irre- 
versibilities in the magnetotransport of mesoscopic multi- 
graphene samples 9 - are compatible with the existence of 
superconducting patches connected by semiconductor re- 
gions, similar to those irreversibilities observed in granu- 
lar high-temperature superconductors 3 ^. Transport mea- 
surements at different regions of the same micrometer 
small and a few tens of nanometers thick multigraphene 
samples reveal an inhomogeneous behavior compatible 
with the existence of metallic- and semiconductor-like 
regions^. Compatible with this view of HOPG, elec- 
tric field force microscopy (EFM) measurements on its 
surface indicate the existence of large variations in the 
electronic potential providing a clear hint that HOPG 
is an inhomogeneous electronic system 3 ^. A recently 
done correlation between the thickness dependence of the 
resistivity and magnetotransport with the internal mi- 
crostructure of HOPG suggests that the internal inter- 
faces between crystalline graphite regions are the regions 
where superconductivity can be located^ 3 -. These inter- 
faces may have enough hydrogen trapped to increase the 
electronic density, triggering superconductivity in differ- 
ent regions with different critical temperatures, keeping 
the quasi-two dimensionality in agreement with the large 
anisotropy observed in the experiments 9 -. 

We should mention other approaches that deal with su- 
perconductivity in graphenc and graphite. Starting with 
graphen o 39 ' 40 we would like to note that for a carrier den- 
sity n < 10 13 cnr 2 the critical temperature T c < 10~ 5 K. 



Therefore, one needs to increase drastically the electron 
concentration in order to have a T c in the few Kelvin 
range. Another publication 41 obtains very large num- 
bers for T c in graphite. The basis of this result is to add 
certain atoms to graphite, in particular sulfur. The at- 
tractive interactions between electrons are based in elec- 
tronic correlations effects. It is basically the same type 
of approach that we use but with phonons. It is clear 
that to produce superconductivity in graphite one has to 
incorporate a material into its planes to reise the carrier 
density to ~ 10 14 cm~ 2 . Otherwise we find that no ap- 
preciable superconductivity can be expected in graphite 
within the mean field approach. 

IV. CONCLUSION 



In conclusion, the superconducting-like behavior ob- 
served in bulk HOPG as well as in mesoscopic multi- 
graphene samples could be explained applying the known 
techniques within the BCS approach and calculating 
A(0) for the quasi-2D problem. From our results we ex- 
pect clear variations of T c with the electronic density, in- 
dicating that neither perfect HOPG nor perfect graphene 
could be superconducting. We have presented results 
that indicate a huge increase in T c for the 2D case with 
respect to the 3D for the same set of parameters. This is 
due to the behavior of QtC^s) and Qf{ts) providing much 
higher values for A(0) in the 2D case. We have applied 
the BCS approach to the system graphene-hydrogen and 
have found that A(0) > 25 K are possible with reasonable 
electronic densities. Upon effective mass and electronic 
density, high critical temperatures in HOPG-hydrogen 
system may be realized and therefore the problem of su- 
perconductivity in graphite should be taken with more 
attention. 
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